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Consider a linear Brownian motion (Bg, s > 0, Bq = 0), and a given drift /x. Exponential functionals 
of the following form 

4^^)= /*dse-2(^« + ^^) (1) 
Jo 

have recently been a subject of common interest for mathematicians and for physicists. 

Some recent mathematical studies have been partly motivated by continuum-time t finance models 
in which most stock price dynamics are assumed to be driven by the exponential of a Brownian motion 
with drift |12| [14|. In such studies, the representation 

e-(^* + /^*) = i?(/J<i«e-2(^« + ^^)) (2) 

where {R{u),u > 0) is a Bessel process, i.e. an element of an important class of diffusions, exhibits 
the importance of the functional a[^^ . Formula (§) is a particular instance of the Lamperti relation 
which expresses (exp(^t),i > 0) as 

e^t =x( f ds e^^"] (3) 



where {X{u)^u > 0) is a semi-stable Markov process (see Q [10| for some applications, partly in 
mathematical finance). 

In physics, these exponential functionals play a central role in the context of one dimensional 
classical diffusion in a random environment. The random variable can indeed be interpreted as a 
trapping time. Its probability distribution controls the anomalous diffusive behaviors of the particle at 
large time in a infinite sample Q |15]. The distribution of A^lf^ occurs when studying the maximum 
reached by the Brownian particle in a drifted Brownian potential |17|. The functional A^jf^ also arises 
in the study of the transport properties of disordered samples of finite length L [23| |24] |21|. 
In fact, A^l^^ represents the continuous space analogue of the random series introduced by Kesten et al. 
[ p!9| ] for the so called " random random walk" . This random series is generated by a linear recurrence 
relation with random coefficients and therefore constitutes a discrete random multiplicative stochastic 
process. a\-^^ also represents a very simple case of a continuous random multiplicative stochastic 
process [ pT| ] which may be related to hyperbolic Brownian motion |29] 

In this article, we discuss some properties of these functionals, concentrating ourselves mostly on 
the mean- value E{\n.A\^'^), in relation with another physical interpretation inspired by the statistical 
mechanics of disordered systems. In these systems, the partition function Z is a functional which 
depends on a set of "quenched" random couplings. In order to obtain the thermodynamic properties 
of the system, one has to compute the average over the disorder of the free-energy F 

E{F) = E{-kT\nZ) (4) 

This calculation can rarely be done exactly. The determination of the probability distribution of F is 
a still more difficult task. One of the very few cases for which such a calculation can be done exactly is 



the Random Energy Model [13|. It is therefore highly desirable to investigate other explicitly solvable 
cases |2^ Q where the usual tools of disordered systems such as replica methods and variational 
techniques can be tested 0. 
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1 Physical motivation : A toy-model for disordered systems 



Let us consider a particle confined on the interval < x < L and submitted to a random force 
F{x) distributed as a Gaussian white noise around some mean value — /q. The corresponding random 
potential is then simply a Brownian motion with drift which can be written in terms of the Wiener 
process as 

U{x) = - r F{y)dy = fox + ^B^ (5) 

JO 

For a given sample, that is for a given realization of the potential U (x), we define the partition function 

Zl= r dxe-l^^i'') (6) 
Jo 

where as usual (3 is the inverse temperature of the system. It is convenient to introduce a = and 
the dimensionless parameter /x = ^ to rewrite as 

Zf = dx e-(«^^ + V^^^) (7) 
Jo 

Therefore, for a = 2, and A^^^ coincide. For a 7^ 2, using the scaling properties of the Brownian 
motion, one obtains 

aL 

(^) (law) 2 f^^^^_2^^^ + B,) (8) 
a Jo 
Hence 

Zi^ a^-) 



a "2" 



The thermal average of any function g{x) of the position of the particle for a given sample will be 
denoted by an upper-bar 

For instance the thermal average and variance of the position read 

More generally, the generating function of the thermal cumulants of the position reads 

In (e^^ = In Z^''^ « ^ - In z[^^ (13) 

These relations show that the statistical properties of the position of the particle in the case = 
require in fact the knowledge of the partition function with an arbitrary drift /x. 
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The fundamental quantity for the statistical mechanics of this disordered system is the free energy 
F^^^ related to the logarithm of the partition function 

i^'^) = -iln4^) (14) 

This paper deals essentially with the statistical properties of the free energy, and particularly with the 
mean free energy over the disorder denoted by 

E{Fi^^)=-^E {in Zi^^) (15) 

This work also gives us the opportunity to bring together various results which were scattered both 
in the physics and mathematics literature. A comparison between them often yields some interesting 
identities with a non trivial probabilistic content. 

The paper is organized as follows. In section ^, we consider the case of random potential with 
zero-drift ^ = 0. We first give the probability distribution of the free energy. We then establish various 
formulae for its mean value, using in particular the Bougerol identity, and show that the replica method 
gives the correct result. We also discuss the asymptotic behavior of the mean free energy in the limit 
of a large sample L — > oo. In section |3|, we discuss the properties of the free energy in the case of a 
random potential with a positive drift fi > 0. We also establish some relations between E (inZ^^^ 

and E (^'^uj^- section ^ we discuss the case where the length of the sample is an independent 
random variable which is exponentially distributed. 

For notations and properties of the special functions appearing in this paper, we refer the reader 



to Lebedev l2Cl. 



2 Case of random potential with zero drift /i = 

2.1 Distribution of the partition function and associated free energy 

The expression of the generating function of Z^^ has already been derived in another context through 
the resummation of the series of moments |23] and by a path-integral approach [^] 
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71" Jo 



oo 

dk cosh — Kik 



(16) 



We also refer the reader to Q, where it is shown that this result may be derived from the Bougerol 
formula (see later 31). To invert the Laplace transform (16), it is convenient to start from the integral 
representation 



a 



oo — 2i/ — cosht 

dt cos kt e V a 







and to perform the integration over k in (|T^) to obtain 

,2 



y/noL 



, irt \ -2a 

dt cos I — ^ j e OiL e 



■ cosh t 



(17) 



(18) 
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We may then use the following identity 



e 



n fP^ , , COsh^t 

-2\/-cosht cosht /-"^ dZ _„z ^ 

' — e e aZ (19) 



a 



Vvra Jo Z2 

to cast ( |l8|) into the form 

.(0) 



(^e-^'^i )=l"dZ e-P^ {Z) (20) 



where 

^2 cosh^ f 



ipV>{Z) = ^— / dtcosht cos — e aL e aZ (21) 

TTQV-^^ Z2 JO \aLJ 

denotes the probability distribution of the partition function z'j^\ A simple change of variables then 
gives the probability distribution (F) of the free energy F^^ = — In Z^^ 

^ B e / cosh^A ^gj. 

p(°)(F) = ^^e2^ r^^^o^ht cosf4)e"^e V « / (22) 
^ TTa^/L Jo \aLj 

It is interesting to compare with a similar formula given in Q . 

In the limit of large L, the probability density ^l(C) of the reduced variable ^ = ^^2!!° *° 
the Gaussian 

1 

This asymptotic result may in fact be directly obtained since from the definition (|^ 

Zf = dx e+^B. ^ L e^B, (24) 

JO JO 



In = + In f e^^^) ^ (25) 



hence 
Since 

'2a sup 







/ pi I \ v2asupi?s /, N 

In ( / dse^^^^) ^ ^ In e ^ ^\N\ (26) 

V JO / L^oo 



where is a normalized Gaussian variable, one has 



1 j^^(o) (1-) (27) 



2.2 An expression of El \nZ^^ ) derived from the generating function 



-P _ e~P^L 



The Frullani identity 

inzf = r^*^ 



(28) 
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may be used to compute the mean of the logarithm of from the generating function of equation 
(|T6[). Using the intermediate regularization 



e[ ln(z5°^) ) = lim 



dp ^ E{ e' 



one obtains 



dk 



2 f°° 
\ / -K JQ hp' 

where C = —T'(l) is Euler's constant. 



r(6)- 

1 - e~'^^^\kcoth{Trk) 



+ C - In a 



(29) 
(30) 



2.3 Bougerol's identity 



Bougerol's identity Q is an identity in law relating two independent Brownian motions {Bs,s > 0) 
and {'~fu,u > 0), and involving the exponential functional we are interested in. The statement is that 
for fixed t 

(law) 



smh.{Bt 



where 



At 



ds = 



(31) 



In the Appendix, we present a simple proof of this identity due to L. Alili, and we refer the reader to 
for further details and possible generalizations for the case of a non-vanishing drift /U 7^ 0. 
Scaling properties of Brownian motion then give 



At 



(law) 



1 

t I du e 




(32) 



It follows that 



sinh(i?t 



(law) 
(law) 



/o^ du e-2\/^^n 



It 



It 



(33) 
(34) 



E\ InA 



E\ In 



Bt 



Starting from the partition function 



X 



Z 



+ \nt 



(0) {^^) 



+°° dx ~Tr: , /sinh(x) 
e in 



'2-Kt 



- A 

a 



aL 



V X 



+ lnt 



(35) 



(36) 



we thus get a new expression for the mean value of the logarithm 



i?(^lnZf 



+°° dx r , / sinhfx) 

e aL In ' 



D VvtqL V X 

Comparison with our previous result (|^) leads to the identity 



+ ln(L) 



(37) 



X 



dx J /sinh(x)\ 2 1'+°° dk 

a/ttoL \ X J tt J -00 k"^ 



1 _ g-^^'^V/t coth(7rA;) 



+ C - In(aL) (38) 
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In order to understand the meaning of this identity, we differentiate both sides with respect to L and 
use the heat equation 



X 

e aL 



d / 1 
dL \ ^JnaL 
After an easy integration by part we obtain 

dx 



1 



X 

e CiL 



d'^ J / sinh(x) 
d'^x \ X 



a 



4 d'^x WvraL 



r+oo 

Q / dk e' 

-oo 



X 

aL 



■aLk 



kcoth.{iTk) 



(39) 



(40) 



This identity is a particular instance of Plancherel's formula, since the Fourier transform of the function 



d"^ ^ / sinh(x) 
d'^x V X 



1 



1 



sinh (x) 



can be easily obtained after an integration in the complex plane 



dx e 



ikx 



1 



1 



sinh (x) 



TT 



k coth ( —k 
2 



(41) 



(42) 



This formula is also encountered as follows in the study of the Hilbert transform of Brownian 
motion 

ds 

If n(de) denotes the characteristic measure of Brownian excursions, and e the generic excursion with 
lifetime y(e) = inf{t > : e{t) = 0}, then the law of Hy = Jq J under n is and formula (|4|) 
may be interpreted as 

TT I — e"'"' n| i - e ^ Hy = x] (44) 



^(Acoth(^) -|A| 



+°° dx ( ^ 

— e n I 1 - e 2 

00 X 



with 



nf 1 — e 2 



Hy = X 



ex 
2 



sinh(f^) 



(45) 



in the particular case 9 = 2. 



2.4 Mean free energy through RepHca method 

Let X = The replica method is based on the identity 

^[X"] - 1 



E[lnX] = lim 



n 



(46) 



In many applications, one proceeds by looking for an analytic continuation n ^ of the expressions 
of integer moments of X. This procedure is in general mathematically ill-founded, but in the present 
case non-integer moments can be computed. 

The integer moments of the partition function are well-known [|l2[ 



fe=i ^ ' 



(47) 
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Since this expression contains a sum of n terms, it has a priori no meaning for non-integer n. However, 
the following integral representation for the moment of order n ^ 

^2 



E 



dx 



r(n + i)7o yJiiaL 



aL 



( sinh^ X \ 



\ a 



(48) 



is valid for any positive real n > 0, as can be shown using the consequence ( p^ of Bougerol's identity 



The following expansions in n, as n ^ 

/ sinh^ x \ " 



1 + nln 



/ sinh^ 



a 



a 



+ o{n) 



r(i) 



1 - n^-^ + o(n) = 1 + n(C + 2 In 2) + o(n) 



r(n + i) 

where C = — r'(l) denotes Euler constant, lead to the expression 



E 



In(zf) 



lim ■ 



E (Z 



{0)\n 



n 



(49) 
(50) 



(51) 



C + 21n2-lna + 



dx r 



e aL In (sinh x) 



(52) 



D V TiaL 

This formula is of course directly related to the previous expression (^) obtained using Bougerol's 
identity since 



X 

+ CO ^™ 

e aL \n{x^)+\n{L) = C + 2\n2-\n a 

oo y/T^aL 



(53) 



2.5 Asymptotic expression of the mean free energy for large L 

The various expressions obtained previously (pC|-37-5^) yield the asymptotic behavior 

1 \ 



E 



In(zf) 



3 



2\l— + C-lna ^= + O , 



as L — > oo 



(54) 



which agrees with 

We now show how the first two terms may in fact be recovered through a direct asymptotic analysis 

of 



Z 



(0) 



dx e 



Using the scaling properties of Brownian motion, one has 



(55) 



(56) 
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It is convenient to set A = 2aL and to define 

Ia = E 



Introducing the one-sided supremum S\ = sup B^, one has 

x<A 



In / dx e-^^ 



(57) 



/a = E{Sa) + E 



In dx e-('^A--Bx) 
Jo 



It follows from the scaling properties of Brownian motion that 



Ia = VAE{Si) + E ln{A f dx e-^('^i " ^ 

Jo 



(58) 



(59) 



Recall that from the reflection principle 5*1 |^i|- This implies 

E{Si) = E{\B^\) -- 

For the remaining part of I a, theorem 1 of Pitman and Yor [^] gives the convergence in law 



Jo A—>oo 



(60) 



(61) 



where Ti and Ti are two independent copies of the first hitting time of 1 by a two-dimensional Bessel 
process starting from 0. Hence 



E 



In (A / dxe-'^i^^-^- 
Jo 



A^oo 



E 



in4(ri + ri^ 



(62) 



The right-hand side of last equation can be evaluated either by direct calculation or by appealing 
to recent results related to the "agreement formula" for Bessel processes [^6|. The direct calculation 
proceeds as follows : one starts again with the identity 

du 



E 



in(ri + ri) 



u 



-u 



and uses the expression of the Laplace transform |18| 
After an intermediate regularization, one obtains 



/o(V2u) 



E 



ln(Ti+fi) 



lim 

e^0+ 



°° du 
u 



lim 



-U 



C - In 2 ^ , , 

2 /o(e) 



Ko{e) 



C- ln2 



(63) 
(64) 

(65) 
(66) 



The other method relies on the following consequence of the "agreement formula" for Bessel processes 
( equation (37) of [p^]) 



Es 



(ri + fi)2- 



1 

2^T(rTM) 



(67) 
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Here 5 = 2(1 + fi) is the dimension of the Bessel process, and Eg denotes the expectation with respect 
to the law of this process, starting at 0. ( For 5 = 2, or equivalently ^ = 0, we simply note E instead 
of £^2)- By differentiating both sides of (|67[) with respect to /i at /i = 0, one recovers (|66[) 



E 



ln(Ti+ri) 



C7-ln2 



Finally, combining (^) and (|68|), we obtain 



Ia 




TT A^oo 



C + ln2 



We therefore recover the first two terms of (^) 



(68) 



(69) 



E 



ln(Z 



(0)^ 



2i/ — + C-lna) 

TT / L^oo 



(70) 



3 Case of random potential with drift /i > 



3.1 Distribution of the partition function Zj^^ and associated free energy 



The probability distribution ip^^\z) of the partition function Z^' has already been obtained in 



(m) 



another context as the solution of the Fokker-Planck equation |21] |11| 



d4^\z) 
dL 



d_ 
dZ 



aZ' 



d4^\z) 
dZ 



+ {{fi + l)aZ-lH[^\z) 



(71) 



satisfying the initial condition 'iP£=q{Z) = 5{Z). The eigenfunction expansion of the solution exhibits 
the following relaxation spectrum with the length L 



^i'\z) = a e 

0<n<| 

aL 



-aLn(/x-n)(zir(/i^ ( ±y^'''\,-2n ( J_ 
r(l + //-n) \aZ) " \aZ 



1 

e aZ 



(72) 



+-—7 / ds e 4 "^ssinhTTs 
47r2 7o 



1 h l- 

2 2 



where are Laguerre's polynomials and Wp^y are Whittaker's functions. Contrary to the case /i = 0, 
there exists for /u > a limit distribution as L ^ 00 

1 



a 



T{^) \aZ 



1 \^+^ -- 



(74) 



In the physics literature, this limit distribution was first obtained in the context of one dimensional 
Brownian diffusion in a Brownian drifted potential, where z'^ plays the role of an effective trapping 
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time |6|, and was recently rediscovered in another context |22|. In the mathematics hterature, this 
hmit distribution 

(75) 



(^) (la^) J_ 
— , — , 



where is a Gamma variable with parameter ^ 



p(Yf, G iy,y + dy) 



2Y„ 



dy 



(76) 



was first obtained by Dufresne |12], then shown in [^] to be another expression of the law of last 
passage times for transient Bessel processes. It is interesting to point out that equation (71) first 
appeared in a paper of Wong |^8| (see also |jl^). However, the precise connection with our work is 
beyond the scope of this paper. 



3.2 Mean free energy 

The Laplace transform of Z^^'^ corresponding to the probability distribution ([7^|73D reads |21] 



aL 



0<n<f 



-aLn(u - n) 2(/x - 2n) f pY^\ 



H 7j [ ds e 4:^^ ^ssinhvrs 



1 h i- 

2 2 



a/ \ M a 



P 



We use again the Frullani identity 



E InZ 



L 



roo 
lO P 

to obtain after an intermediate regularisation as in (|2 



E InZ 



in a — — y — e 

Flu ^-^ „ nia — n) 



^—aLn{fi — n) 

n) 
aL 



sinhvrs -— r(/^^ + ^^^ 

e 4 



ds 



10 fi^ + s'^ cosh vrs — cos vr/i 
Obviously, the two first terms correspond to the contribution of the limit distribution 

r'(^) 



E[ InZ^Ji^ = -Ina 



Now, it is easily deduced from (|3.2D and (p^) that 



ii;(^lnZ(^)j-i?(^lnZi'')j^^^<! 



r(/^) 



—aL 



(77) 
(78) 

(79) 

(80) 
(81) 

(82) 



1 

^/ttoL 



if ;u > 2 
if ;u = 2 



(83) 



/i^(l — cos TT/i) \aL 



3 a/U r 

vr \ 2 ; — L 

e 4 if ^ < 2 
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It is interesting to compare these results to those obtained in |17| for the hmiting distribution of the 
maximum of a diffusion process in a Brownian drifted environment. The latter may be stated as 

>-2A „-2(/.-l)L 



E 



K 



e 
-2L 



L 



li>2 
\i IJi = 2 



(84) 



where A^^^ and are two independant functionals of type A^^^ defined in (|^). The constant K 

given in [17| [| as the following 4- fold integral 
2I-M 1 



K 



TT r(^) Jo 



dy / dz e-i^^+y '^ 







da a 



z 







z + a 



e 2 







may in fact be explicitely computed to give in the end the simple result 

\2 



1 



nsinhn 6-^2'"°''^" (85) 



(86) 



r(;u) (1 — cos vr^) 

For ^ > 2, the result ( p^ ) therefore coincide with the result (^), where for the particular value a = 2 
Z^^^ reduces to J^^^ (see (^) and (|l|)). To understand this coincidence, we write 



where is a variable distributed as and independent of J^l^ . Therefore 



(87) 



e{ \^A^m -e( InA^'^A =e( In 



1 + e" 



-2(Bi + //L)A^^ 



A 



(m) 



£; In 



1 + 



-^00 



A 



The comparison between (pH) and (pi) therefore leads to 



E(^ln 



-^00 



A 



~ E 



^00 



for fi>2 



This is likely to be understood as a consequence of the following plausible statement 



If 



Xn 



(a.s.) 



0, then E 



ln(l+X„ 



E 



(89) 



(90) 



l + X„ 

which presumably holds for a large class of random variables but the precise conditions of 

validity of ( pOD elude us. However, (^) does not hold for /i < 2 since, in this case, the prefactors in 
(H) and (g) differ. 

To go further into the comparison, we have computed exactly the quantity occuring in ( p^ for 
arbitrary L. We start from the identity 



E 



■^00 



dp E e-P^^^^ ) E (^e-P^L 



(91) 



^In fact (|8^) corrects the formula for the constant K found in The need to divide by ^/n is due to a misprint in 
j], where on page 528, just after formula (6.e), one should read 



(27r3t)l/2 



■ ■ instead of ■ 



(27r2t)l 
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Using the following consequences of ( [78|) 



and 



aL 



M-l 
P \ 2 



^M-i(2\/- 



we get 
E 



yip) 



2 2 



2 , 2^ 



„ / dxKa-iix) / ds e 4 ssinhvrs 



2 2 



For /i < 2, the order of integrations may be exchanged to give 

aL 



E 



1 



which reduces for ^ = 2 to 
E 



2r(/x) Jo 



7(2) 



ds e 4 



s sinh TTS 



cosh TTS — cos TT// 



i h 

2 2 



(92) 

(93) 

Kis{x) 
(94) 

(95) 



7(2) ^(-2) 



2 Jo 



aL 
dse 4 



(4 + s^) vrs^cosh^ 



sinh^ ^ 



(96) 



For fi > 2, one cannot exchange the order of integrations in (p^). However, one may start from a 
series representation ( eq (5.9) in pl|) of the generating function ([9^) to obtain, after some algebra 



involving deformation of a contour integral in the complex plane (see |21] for a similar approach) the 
general result for arbitrary /i > 

"r(^) 



E 



J2 (^-2n)^^^^^^^^~"^ ,-aLn{fi-n) 



(97) 



l<n<\ 

aL 



+ 



2r(//) Jo 



ds e 4 



(^^ + s^) ssinh 



TTS 



cosh vrs — cos vr/i 



i h 

2 2 



(98) 



From (05-96-08), one easily recovers the corresponding asymptotic results of (p4), which were obtained 



in [17| by a quite different method, relying on the computations made in [^]. The presence of discrete 
terms for ii> 2 again explains the transition at = 2 of the asymptotic behavior. 



3.3 Some relations between iJ^lnZl^^^ and mean inverse e(^-^ 
The first negative moment can be obtained from the generating function (|7^) written above 



a ^ (// — 2n) e 



-aLn{n — n) 



+ 



a r 
2 Jo 



ds 



s sinh vrs 



aL 



cosh vrs — cosTT/u 



1) 
(99) 

(100) 
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The previous explicit expressions therefore lead to the very simple identity for any /i > 

A£(,„4-)).£(_L)_„, (101) 

Can it be derived directly using only basic properties of Brownian motion ? Trying to do so gives 
in fact two other identities of the same kind, but not ( [L01[ ). The first one relates the exponential 
functional for two opposite drifts and (— ^) 

f,E{>nZf)=E{^) (102) 

This identity can be obtained from a simple reparametrisation x' = L — x in the denominator of the 
left handside of expression 



^^f 1„ .-("^- + ■^B.A = E ^-<°''"+^^°^' ] (103) 

dL \ Jo J ^-{afxx + V2aB^) J 

The second one relates the exponential functionals for two dimensionless drifts differing by two 

|-ij(l„Z<''))=e-"i(*'-l)£(-J^) (104) 

L 

This identity follows from Cameron-Martin or Girsanov relations. In fact these relations lead to a 
more general relation between the two characteristic functions of the exponential functionals for (/x) 
and — 2) 

^E^e^P^''^^ = -p e-«^(/^ - 1) E(e-<'~''^ (105) 



4 Expressions of E'^ln [Zj^^^f) with an independent exponential length 

Lx 



It is well known that the laws of additive functionals of Brownian motion with drift /j, 

A{ ^=^^ f'dsfiBs+ns) (106) 
Jo 

may be easier to compute when the fixed time t is replaced by an independent exponential time T\ of 
parameter A 

p( Txe[t,t + dt]) = Xe-^^dt (107) 



This is indeed the case for A{ = a[^^ i.e. f{x) = It was shown in pO] p3|| that 



^'^=^ ^ where ^ = V2AT^ V2A + + /. 
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^a,i3 denotes a Beta- variable with parameters (a, /?) 



P[ Xo,pe[x,x + dx\ ] = — dx (0<x<l) (109) 



and denotes a Gamma- variable of parameter (7) (iTq). Consequently, one has 

Vv ; 2"r(l + a + n)r(6) ^ ^ 

and 

- £ 



E( InA^^^ ) = C7 + ln2 + ^(1 + a) +^^(6) (HI) 



where C = — r'(l) is Euler's constant, and ij{x) = yJxJ- Classical integral representation of the 
function ip allows to invert the Laplace transform in A implicit in ( |111| ), hence to recover e( In^'^^V 



however the formulae we have obtained in this way are not simple. 

To transpose the result (111) for the partition function Zj^^ describing the case where the length 
of the disordered sample is exponentially distributed, one only needs to use the identity derived from 
(I) 

(^) (law) 2 (^) 



Appendix : A simple proof of Bougerol's identity 

As Bougerol's identity (^) may appear quite mysterious at first sight, we find useful to reproduce 
here a simple proof of this identity due to L. Alili and D. Dufresne. We refer the reader to and [||] 
for further details and possible generalizations for the case of a non-vanishing drift ^ 7^ 0. 

Consider the Markov process 

Xt = e^* / e-^= djs (113) 
Jo 

where Bt and 7^ are two independent Brownian motions. Ito formula yields the stochastic differential 
equation 

dXt = ^Xt dt + {Xt dBt + djt) (114) 
We introduce a new Brownian motion /?( by setting 



Xt dBt + d^t = \JX^ + 1 dPt (115) 

from which it follows that 

I 1 
dXt = -Xt dt + [Xj + 1) ^ d(5t (116) 

The comparison with 

d[sinh(A)] = ]-[sinh{(it)]dt + {sinh\(5t) + l) ^ d(5t (117) 
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shows that 

sinh(A) ^ Xt = e^* t e-^= ^7. (118) 

Jo 

The use of scaUng properties of Brownian motion finally gives 

sinh(A) ^'=^ 7,(0) with A^"^ = f e^^' ds (119) 

Jo 

where 7 denotes a Brownian motion, which is independent of Bs- 
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